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this  reason  a  theory  of  filtration  has  been  developed  and  is  presented 
herein.  This  is  the  only  theoretical  treatment  applicable  to  tribology,  al¬ 
though  an  alternative  purely  empirical  model  due  to  Fitch  is  also  available. 

In  contrast  to  the  empirical  approach,  the  results  are  not  dependent  on  a  spe¬ 
cific  particle  size  distribution  ("AC  Fine  Test  Dust”),  and  the  effect  of 
governing  parameters  (filter  grain  size  and  void  fraction,  fluid  flow  rate  and 
viscosity,  etc. )  can  be  predicted.  The  predictions  of  the  theory  are  consist¬ 
ent  with  physical  reasoning  and  the  known  properties  of  filters. 

The  basic  concepts  developed  can  be  directly  applied  to  Naval  lubrication 
systems,  with  a  very,  quick  pay-back  in  terms  of  reduced  maintenance  and  reli- 
ability  costs. it  is  shown  that  the  optimum  filter  is  an  "absolute11  filter 
which  traps  all  particles  greater  than  some  size  aQ  and  allows  all  smaller  ones 
to  pass  through.  Since  real  filters  are  not  absolute  but  possess  a  character¬ 
istic  efficiency  curve,  the  theory  can  be  used  to  select  the  filtration 
parameter  values  which  will  result  in  the  most  nearly  ideal  performance.  The 
ideal  "cut-off"  size  Sq  can  be  determined  from  a  simple  series  of  field  tests 
which  could  be  readily  implemented  as  part  of  Naval  ship  procedure. 

Experimental  research  to  date  indicates  that  the  particle  size  distribu¬ 
tions  of  actual  contaminated  lubricants  vary  greatly  and  cannot  be  character¬ 
ized  solely  as  AC  Fine  Test  Oust.  It  appears  that  filterable  size  particles 
(>  5  Jim)  in  used  oil  have  little  effect  in  the  rheological  flow  properties  of 
the  oil.  The  well-known  viscosity  level  changes  are  probably  due  to  chemical 
contamination  and  sub-micron  soot  particles.  These  findings  may  have  a  bearing 
on  efforts  to  increase  the  useful  life  of  circulating  lubricants. 

Future  research  efforts  are  suggested.  It  is  proposed  that  the  filtration 
theory  be  continued  to  include  additional  effects,  and  to  put  forth  the  results 
in  a  user-accessible  program.  In  addition,  filtration  experiments  will  be  con¬ 
ducted  to  check  the  theory  in  both  a  highly  controlled  laboratory  situation, 
and  in  a  more  real-life  setting.  Finally,  it  is  proposed  that  a  program  of 
'lubrication  system  filtration  optimization  tests  be  conducted  as  suggested  in 
the  text.^  Jfie  goal  here*  is  to  set  up  procedures  such  that  optimum  filtration 
can  be  achieved  in  Naval  ship  lubrication  systems.  _ 
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NOMENCLATURE 


a(-  “g”")  "  Viscous  Effect  Number,  Eq.  (16) 

a  -  Particle  radius 
a^«  -  Relative  particle  radius 

aQ  -  Absolute  filtration  cut-off  size  -  Pigs.l  and  14,  Eq. (38) 

B, C,D  -  Auxiliary  parameters,  Eq. (A. 4) 

d  (*  2a)  -  particle  diameter  or  major  dimension 

— i 

Fr,Fg  -  Viscous  force  (in  the  r, 8-direction)  due  to  particle 
motion.  Superscript  notation  in  the  Appendix,  p.A2. 

Bars  denote  dimensionless  values,  see  Eq. (A. 7) 

f*,f^  -  Viscous  force  (in  the  r, 9-direction)  due  to  fluid  motion. 

Superscript  notation  in  the  Appendix,  p.A2.  Bars  denote 
dimensionless  values,  see  Eqs. (A. 7) 

~v,Fvr,Pv0  ”  Tot®!  force  on  a  particle  due  to  viscous  fluid  effects? 

~  vector  form,  r, 9  component 

P1,Plr,P1g  -  Force  on  a  particle  due  to  particle  translation;  vector 
form,  r, 9-component 

-  Force  on  a  particle  due  to  particle  rotation;  vector 

~  form,  r, 9-component 

~1J ^lr' f19  ”  F°rce  °n  a  particle  due  to  fluid  shearing;  vector  form, 

~  r, 9-component 

f2, fjr*  f29  ”  Force  on  a  particle  due  to  fluid  vorticity;  vector  form, 
r, 8-component 

gl,g2,g3  “  Auxiliary  parameters,  Appendix,  Eq. (A. 10) 

13 

H  -  Hamaker's  Constant  (~  10  ergs).  Determines  magnitude 
of  van-der-Waals  force,  see  Eq.  (A.  7) 

K^,K2,K3,K4  -  Auxiliary  parameters,  Eq. (A. 2) 

L  -  Filter  length 

m  -  Particle  mass 

n,n^  c  ^  -  (Total)  number  of  particles  per  unit  volume.  Subscripts 

*  *  *  denote  total,  captured,  inlet  and  exit,  respectively. 

/  L  \ 

N  »  )  “  F!!f®r  length  number,  number  of  unit  cell  layers 


p  -  Auxiliary  parameter,  Eq. (A. 1) 


r.r 


Radial  coordinate,  radial  coordinate  of  particle. 
Bus  denote  dimensionless  values,  see  Eq.  (16),  Fig.  3 


R  -  Radius  of  collector 


R  -  Radius  of  unit  cell 
s 


2 

9  |1R  / 

T1 

9 


Reynolds  number 


Stokes  number 

Viscous  torque  (in  the  cp-direction)  due  to  particle  motion. 
Superscript  notation  in  the  Appendix,  p.A2.  Bars  denote 
dimensionless  values,  see  Eqs. (A. 7) 


t1  -  Viscous  torque  (in  the  cp-direction)  due  to  fluid  motion. 
^  Superscript  notation  in  the  Appendix,  p.A2.  Bars  denote 
dimensionless  values,  see  Eqs. (A.7). 


t,t 


e,z,x 


Time,  dimensionless  time,  Eq.  (17) 

Particle  velocity  component  in  the  r, 0,z  or  x-direction. 
Bars  denote  dimensionless  values,  see  Eqs. (16)  and  (A. 2) 


V  -  Velocity  of  fluid  entering  filter 


0,z,x 


Fluid  velocity  component  in  the  r, 0,z  or  x-direction.  Bus 
denote  dimensionless  values,  see  Eqs.  (16)  and  (A. 2) 


W  -  Auxiliary  puameter,  Eq.  (A.l) 


51, 5,  z 
’  *  P 

P 


e 


<D,U) 


Alternative  dimensionless  coordinates,  Eq. (A. 2) 
Filtration  beta  ratio,  Eq. (2) 

Porosity,  volume  void  fraction 

Angulu  coordinate,  angular  coordinate  of  puticle 
Density  of  fluid,  puticle 
Fluid  viscosity 

Puticle  rotation  rate,  dimensionless  puticle  rotation 
rate,  Eq. (16) 
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ABSTRACT 


Proper  lubrication  filtration  is  essential  if  maintenance  and  reliability 
costs  of  mechanical  systems  are  to  be  minimized.  Filtering  too  little  can  cause 
high  rates  of  wear  and  costly  failures,  while  filtering  too  much  can  cause  ex¬ 
cessive  maintenance  costs.  Better  understanding  is  required  such  that  the  op¬ 
timum  level  of  filtration,  between  these  two  extremes,  can  be  achieved.  Opti¬ 
mum  filtration  is  equally  desirable  even  in  the  presence  of  an  on-line  wear 
monitoring  system. 

For  this  reason  a  theory  of  filtration  has  been  developed  and  is  presented 
herein.  This  is  the  only  theoretical  treatment  applicable  to  tribology,  al¬ 
though  an  alternative  purely  empirical  model  due  to  Fitch  is  also  available. 

In  contrast  to  the  empirical  approach,  the  results  are  not  dependent  on  a  spe¬ 
cific  particle  size  distribution  ("AC  Fine  Test  Dust"),  and  the  effect  of 
governing  parameters  (filter  grain  size  and  void  fraction,  fluid  flow  rate  and 
viscosity,  etc. )  can  be  predicted.  'Hie  predictions  of  the  theory  are  consistent 
with  physical  reasoning  and  the  known  properties  of  filters. 

The  basic  concepts  developed  can  be  directly  applied  to  Naval  lubrication 
systems,  with  a  very  quick  pay-back  in  terms  of  reduced  maintenance  and  reli¬ 
ability  costs.  It  is  shown  that  the  optimum  filter  is  an  "absolute"  filter 
which  traps  all  particles  greater  than  some  size  aQ  and  allows  all  smaller  ones 
to  pass  through.  Since  real  filters  are  not  absolute  but  possess  a  character¬ 
istic  efficiency  curve,  the  theory  earn  be  used  to  select  the  filtration  parameter 
values  Which  will  result  in  the  most  nearly  ideal  performance.  The  ideal  "cut¬ 
off  size  aQ  can  be  determined  from  a  simple  series  of  field  tests  which  could 
be  readily  implemented  as  part  of  Naval  ship  procedure. 

Experimental  research  to  date  indicates  that  the  particle  size  distribu¬ 
tions  of  actual  contaminated  lubricants  vary  greatly  and  cannot  be  characterized 
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solely  as  AC  Fine  Test  Dust.  It  appears  that  filterable  size  particles  (>  5  urn) 
in  used  oil  have  little  effect  in  the  rheological  flow  properties  of  the  oil. 
The  well-known  viscosity  level  changes  are  probably  due  to  chemical  contamina¬ 
tion  and  sub-micron  soot  particles.  These  findings  may  have  a  bearing  on  ef¬ 
forts  to  increase  the  useful  life  of  circulating  lubricants. 

Future  research  efforts  are  suggested.  It  is  proposed  that  the  filtra¬ 
tion  theory  be  continued  to  include  additional  effects,  and  to  put  forth  the 
results  in  a  user-accessible  program.  In  addition,  filtration  experiments  will 
be  conducted  to  check  the  theory  in  both  a  highly  controlled  laboratory  situa¬ 
tion,  and  in  a  more  real-life  setting.  Finally,  it  is  proposed  that  a  program 
of  lubrication  system  filtration  optimization  tests  be  conducted  as  suggested 
in  the  text.  The  goal  here  is  to  set  up  procedures  such  that  optimum  filtra¬ 
tion  can  be  achieved  in  Naval  ship  lubrication  systems. 
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I.  INTRODUCTION 


On  I  December  1979,  a  research  program  at  Rensselaer  Polytechnic  Insti¬ 
tute  by  the  Office  of  Naval  Research  (ONR)  was  begun  entitled  "Control  of 
Contaminants  in  Circulating  Engine  Lubricating  Oil,"  ONR  N00014-79-C-0100, 

Dr.  John  A.  Tichy,  Principal  Investigator.  The  original  proposal  was  struc¬ 
tured  around  a  two-year  program,  although  only  the  first  year  was  initially 
funded. 

Although  the  original  research  program  called  for  four  essential  task 
areas,  by  mutual  agreement  of  Dr.  Tichy  and  ONR  contract  monitor,  Commander 
H.P.  Martin,  the  one-year  program  would  stress  only  the  first  two  tasks, 
these  areas  are: 

(1)  Determine  the  Particulate  Content  of  used  Oils  and 
Correlate  with  Measured  Rheological  Behavior 

(2)  Develop  a  Lubrication  Filtration  Theory. 

this  report  describes  the  research  effort  to  date  in  these  two  task 
areas,  the  principal  investigator  believes  that  the  filtration  theory 
development  has  been  very  successful  and  holds  great  promise  for  a  quick 
pay-back  to  the  Navy.  For  this  reason,  the  bulk  of  this  report  stresses 
this  area,  in  particular,  Section  IV  outlines  the  results  of  the  theory 
itself  and  Section  V  explains  an  application  which  could  be  put  to  use 
immediately  after  a  fairly  simple  on-board  test  program. 

Section  VI  describes  the  experimental  work  to  date  associated  with 
Task  1.  Future  research  recommendations  are  presented  in  Section  VII.  The 
recommendations  deviate  somewhat  from  the  tasks  of  the  original  proposal  in 
light  of  the  research  performed  to  date. 
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II.  THE  NEED  FOR  A  THEORY  OF  LUBRICATION  FILTRATION 


1.  Seeking  the  Optimum  Level  of  Filtration 

1.1  Maintenance,  Reliability  and  Optimum  Filtration 

Recent  studies  of  maintenance  procedures  in  naval  diesel-powered 
ships  tH  and  aircraft  engines  [2]  have  indicated  that  the  contamination  of 
circulating  lubricating  oil  is  a  source  of  increased  wear  and  friction, 
failure  of  machine  elements,  and  high  maintenance  costs.  In  many  areas  the 
costs  of  maintenance  are  equivalent  to,  or  may  even  exceed,  capital  costs. 
Consequently,  maintenance  costs  are  expected  to  play  an  increasing  role  in 
managerial  decision-making. 

The  filtration  of  particulate  matter  from  lubricating  oil  is  an 
important  component  of  the  maintenance  of  any  engine  power  plant  mechanical 
system.  In  these  days  of  advanced  technology,  filtration  by  porous  media 
remains  by  far  the  most  practical  means  to  remove  harmful  contaminants*  from 
an  engine's  lubricating  system,  although  the  basic  process  is  still  poorly 
understood. 

The  filtration  process  may  serve  to  either  increase  or  decrease 
maintenance  costs.  Decreased  costs  may  occur  because  the  reduced  level  of 
abrasive  material  in  the  lubricant  causes  a  reduced  wear  rate.  The  exact 
relation  between  particle  size  and  abrasive  wear  rate  has  been  studied  re¬ 
cently,  at  least  under  laboratory  conditions  (e.g. ,  [3,4, 5,6]).  These  workers 


Contamination  in  lubricants  may  refer  to  the  harmful  effects  of  relatively 
small  particles  (say  d  <  1  pm)  such  as  those  from  dispersed  diesel  engine 
"soot",  or  effects  of  larger  particles  (say  d  >  3  pm)  from  wear  debris  or  the 
ingression  of  dust  and  dirt.  In  the  former  case,  deleterious  effects  are  due 
to  such  things  as  the  build-up  of  sludge,  viscosity  increase,  chemical  or 
additive  property  change,  etc.  In  the  latter  case  the  larger  particles  di¬ 
rectly  cause  abrasive  wear.  Unless  a  prior  distinction  is  made,  the  term 
contamination  refers  here  to  the  latter  condition. 


3 


generally  conclude  that  wear  volume  rapidly  increases  with  the  abrasive  grit 
size  for  15  (Jbm  <  d  <  100  tun  (where  d  is  the  particle  diameter,  or  major 
dimension)  and  then  increases  at  a  much  slower  rate,  or  not  at  all,  for  larger 
particles.  Thus  reduced  wear  rate  due  to  filtration  helps  to  extend  component 
life  and  reduce  costs  due  to  sudden  unanticipated  breakdowns.  On  the  other 
hand,  higher  maintenance  costs  are  associated  with  increasingly  finer  filtra¬ 
tion.  Such  filters  may  be  expensive  themselves,  clog  quickly,  and  require 
costly  down-time  for  replacement.  Apparently,  then,  an  optimum  (economic) 
level  of  filtration  must  exist  for  any  mechanical  lubrication  system,  as  a 
trade-off  between  these  two  conditions. 


1.2  Filtration  Performance  Curves 

The  situation  is  made  more  complicated  by  the  fact  that  filters  are 

not  "absolute,"  i.e.,  they  do  not  capture  all  particles  above  some  size  and 

allow  all  others  to  pus.  Father  each  filter  has  a  characteristic  curve  of 

efficiency  versus  particle  size,  see  Figure  1.  At  each  particle  radius 
* 

size  a  ,  an  efficiency  7]  is  defined  by 


11(a) 


n  (a)  1-n  (a) 

c  e 


n^  (a)  ni(a) 


(1) 


where  n  (a)  is  the  number  of  particles  (of  size  a,  per  unit  volume)  captured 
c 

by  the  filter,  n^(a)  is  the  number  of  particles  at  the  filter  inlet  and  nfi(a) 
is  the  number  at  the  filter  exit. 

How  then  does  the  designer  decide,  for  example,  whether  one  filter 
which  removes  10  iim  particles  at  90%  efficiency  is  justified  over  another  which 


* 

Following  traditional  usage,  it  is  most  common  to  use  radial  size  a  in 
theoretical  work,  and  diametrical  size  (major  dimension)  d  ■  2a  in  applications. 


Efficiency  i? 


Figure  1  Filtration  Efficiency  Curves 
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removes  20  pm  particles  at  60%  efficiency?  It  will  be  shown  that  the  earlier 
statement  can  be  taken  further.  Not  only  does  a  general  "optimum  level  of 
filtration"  exist,  but  a  specific  filtration  efficiency  curve  can  be  shown  to 
be  the  optimum  for  any  given  system. 

1.3  Filtration  Theory  Versus  Filtration  Model 

In  this  report  the  term  "filtration  theory"  is  distinguished  from 
"filtration  model."  The  latter  are  based  solely  on  the  empirically  measured 
behavior  of  filters  under  a  set  of  conditions  -  with  the  results  correlated  in 
some  systematic  way.  Predictive  ability  by  the  model  is  not  based  on  physical 
reasoning  but  only  on  the  notion  that  filters  behave  pretty  much  alike  and  the 
behavior  of  any  new  filter  can  be  expected  to  be  the  same  as  others.  A  model 
of  this  sort  cannot  be  expected  to  predict  behavior  if  actual  performance 
conditions  are  significantly  different  from  test  conditions.  A  filtration 
theory,  on  the  other  hand,  is  based  on  a  physical  concept  of  what  happens  in¬ 
side  the  filter,  i.e.,  how  it  works.  This  physical  picture  is  then  set  in 
mathematical  terms  and  overall  behavior  of  the  filter  is  predicted.  The  effect 
of  different  variables  on  filter  performance  can  then  be  determined.  A  theory 
of  this  sort  has  predictive  ability  under  a  wide  range  of  conditions  provided 
the  original  physical  scenario  is  not  violated.  Due  to  the  large  number  of 
variables,  and  the  complexities  involved,  it  is  unlikely  that  a  completely 
accurate  theory  can  be  currently  formulated.  On  the  other  hand,  even  a  simp¬ 
lified  theory  cm  yield  information  which  can  be  used  as  a  guideline  in  develop¬ 
ing  accurate  correlations  frcn  experimental  data.  Ultimately,  of  course,  any 
physical  theory  must  be  tested  against  real-world  behavior. 
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1.4  Over-Filtration  and  On-Line  Monitoring 

The  current  trend  in  the  Navy,  apparently,  is  away  from  increasingly 
finer  filtration  but  to  better  "tune"  the  filter  to  the  system's  mechanical 
and  economic  needs.  This  less  conservative  filter  design  requires  a  higher 
level  of  knowledge  of  filtration  behavior,  i.e.,  a  theory  of  lubricant  filtra¬ 
tion.  In  the  same  way  less  conservative  aircraft  design  (reduced  safety  factors, 
etc. )  requires  ever  more  sophisticated  analytical  techniques  to  predict  behavior 
more  accuratelv. 

The  Navy  is  also  undertaking  a  program  in  on-line  monitoring  of  the 
lubricating  oil  of  power  plant  systems.  It  is  hoped  that  less  dependence  on 
a  high  level  of  filtration  will  be  required  if  the  machine' s  "health"  is  con¬ 
stantly  monitored.  This  study  can  be  seen  as  a  complement  to  such  a  program. 
Even  if  the  onset  of  machine  element  failure  can  be  determined  accurately  and 
inexpensively  by  on-line  monitoring  there  still  will  be  major  costs  associated 
with  down-time  and  replacement.  Optimal  filtration  will  still  be  every  bit  as 
desirable  a  goal.  In  addition,  there  are  probably  many  small  ship  systems 
where  a  sophisticated  monitoring  program  cannot  be  cost  effective,  but  where 
the  proper  filtration  level  can  significantly  reduce  overall  costs. 

2.  Theories  and  Models  of  Filtration 

2.1  Theories  from  Chemical  and  Environmental  Engineering 

It  is  not  well  known  among  tribologists  that  a  theory  of  filtration 
has  been  developed  over  some  50  years  and  is  currently  used  with  some  success 
in  the  engineering  of  depth  filtration  of  water  and  fiber  filtration  of 
aerosols.  The  theories  are  based  on  porous  media  models  for  the  filter  it¬ 
self,  e.g. ,  an  array  of  spherical  or  cylindrical  collectors.  The  fluid  and 
particles  are  forced  to  follow  a  tortuous  path  through  the  filter  medium  and 
the  filtering  action  is  due  to  capture  of  the  contaminant  particle  when  it 
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collides  with  the  filter  grains  or  fibers  (collectors),  see  Figure  2.  The  total 
efficiency  is  determined  from  a  single  fiber,  or  single  grain  efficiency;  the 
so-called  unit  cell  efficiency. 

Hie  unit  cell  is  generally  modeled  as  a  fluid  flow  past  a  cylindrical 
or  spherical  obstacle  with  an  entrained  particle,  see  Figure  3.  In  the  more 
sophisticated  theories,  forces  and  moments  are  balanced  on  the  contaminant 
particle  and  solution  of  the  resulting  equations  yields  the  particle  trajectories. 
Some  particles  collide  with  the  collector  while  others  escape  to  a  greater  filter 
depth  and  encounter  another  unit  cell.  The  efficiency  is  then  the  fraction  of 
particles  entering  the  filter  which  ultimately  collides  with  a  collector. 

The  basic  theory  described  is  outlined  in  books  by  Davies  [7]  and 
Theodore  and  Buonicore  [8] ,  and  in  a  series  of  articles  by  Tien,  Rajagopolan 
et  al.  [9,10,11].  Many  other  references  are  also  available  but  it  serves  no 
purpose  to  prepare  an  extensive  bibliography  here. 

2.2  Filtration  in  Tribological  Applications 

Other  than  that  proposed  here,  there  is  no  theory  of  lubrication 
filtration.  The  role  of  particulates  and  contamination  on  lubricant  perform¬ 
ance  has  just  recently  been  appreciated.  Only  the  most  basic  scientific 
studies  of  filtration  on  engine  wear  and  bench  test  studies  of  filtration  on 
bearing  wear  and  performance  have  been  performed  to  date  (e.g. ,  [12,13,14,15]). 
Perhaps  this  is  why  attempts  at  rigorous  engineering  study  of  the  filtration 
process  itself  have  been  slow  to  develop,  combined  with  skepticism  that  a 
theory  could  do  much  good  on  so  complex  a  problem. 

In  the  theory  proposed  here  (details  are  presented  in  the  following 
Section  III  and  Appendix  A)  the  basic  concepts  described  in  Section  II,  2.1, 
have  been  adapted  to  tribology  conditions.  Although  the  fundamental  physical 
model  remains  the  same,  the  particular  forces  and  moments  which  are  significant 
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Figure  3  The  Unit  Cell 
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depend  on  the  specific  type  of  application,  as  do  the  range  of  variables 
involved.  For  example,  in  tribology,  viscous  fluid  forces  predominate  and 
particle  sizes  are  in  the  range  (say)  1  pm  <  d  <  100  pm.  As  a  contrary 
example,  in  clean  water  filtration,  contaminant  particles  may  be  0.1  pm  and 
Brownian  motion  diffusion  forces  become  important.  Other  aspects  not  gener¬ 
ally  considered  in  the  chemical  engineering  theories  are  the  effect  of  a 
build-up  of  deposit  (clogging)  and  the  effect  of  particle  inertia.  In  lubri¬ 
cation  filtration,  the  fluid  is  often  vigorously  pumped  through  the  filter 
medium  and  the  usual  creeping  or  Stokes  flow  assumption  (the  neglecting  of 
inertia)  may  be  violated. 


3.  Hie  "Beta  Ten"  Filtration  Model 
3.1  Description  of  the  Model 

Fitch  and  co-workers  [16]  have  developed  a  filtration  model  applicable 
to  systems  where  fluid  is  circulated  and  particulate  contamination  of  this 
fluid  tends  to  produce  harmful  wear.  Fitch' s  work  to  date  is  the  only  pub¬ 
lished  material  which  deals  with  lubricant  filtration  on  any  sort  of  scien¬ 
tific  basis.  Fitch  has  established  a  so-called  filtration  beta  ratio  [17]  as 
the  basis  of  evaluation  of  a  filter's  performance.  Hie  beta  ratio  at  some 
particle  size  a  is  defined  by  the  following  equation: 


n.  (a  >  a  ) 
^(ao5  n  (a  >  a  ) 


(2) 


i.e.,  the  beta  ratio  is  the  ratio  of  the  number  of  particles  (per  unit  volume) 

greater  than  size  aQ  at  the  inlet  n^  (a  >  aQ)  to  the  number  of  particles 

greater  than  size  a  at  the  exit  n  (a  >  a  ).  thus  the  higher  the  B,  the  better 
o  e  o 

the  filtration.  Hie  beta  values  vary  with  particle  size  producing  a  different 
sort  of  performance  curve  than  the  efficiency  curve  described  above. 


T 

i 


Hie  Beta-ten  model  [18]  is  based  solely  on  the  measured  behavior  of 


filters  with  a  prescribed  ingress ion  rate  of  a  particular  contaminant  particle 
size  distribution,  namely,  "AC  Fine  Test  Dust".  Indeed  the  entire  beta  ratio 
concept  and  definition,  Eq. (2),  is  based  on  the  AC  Fine  Test  Dust  particle  size 
distribution.  Beta  curves  (0  vs.  particle  size  a)  of  numerous  filters  were 
plotted  on  special  semi-log  paper.  Each  curve  is  then  specified  by  its  0^Q- 
value  -  the  beta  ratio  for  10  lun  particles.  The  Beta- ten  model  performance 
curve  is  shown  later  in  Figure  13. 

3.2  Shortcomings  of  the  Model 

Perhaps  the  most  notable  contribution  of  the  Fitch  approach  has  been 
to  convince  tribologists  that  filters  are  not  "absolute",  as  described  in 
Section  II,  1.2. 

Although  Fitch*  s  work  has  been  instrumental  in  elevating  filtration 
in  tribology  from  a  black  art,  there  are  nevertheless  numerous  shortcomings 
in  his  approach.  They  are  as  follows: 

(a)  Hie  performance  curve  of  any  filter  with  the  Fitch  model  depends 
only  on  one  parameter,  its  01Q-rating.  Both  the  theory  to  be 
presented,  and  conanon  sense,  would  indicate  that  many  parameters 
have  influence  -  e.g.,  filter  collector  size,  porosity,  material, 
fluid  flow  rate,  viscosity,  length,  particle  concentration,  etc. 

(b)  The  entire  concept  i3  based  on  the  particle  size  distribution 
of  AC  Fine  Test  Dust,  if  a  different  distribution  function 
occurs  in  the  contaminated  lubricant,  and  research  presented 

in  Section  V  suggests  this  is  often  the  case,  the  significance 
of  the  beta  ratio  is  unclear. 

(c)  Effects  of  variable  changes  cannot  be  predicted.  The  model  will 
not  predict  what  can  be  expected  if,  for  example,  fluid  flow 
rate  is  increased,  temperature  is  increased,  filter  length  is 
decreased,  the  filter  is  packed  more  densely,  etc. 
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(d)  Time-dependent  effects  such  as  those  due  to  clogging  cannot 
be  accommodated. 

(e)  The  Beta-ten  model  cannot  suggest  changes  in  the  filter,  or  in 
in  filtering  procedure,  to  obtain  an  efficiency  curve  better 
suited  to  the  lubricating  system. 

It  seems  clear  that  the  work  of  Fitch  and  Tessman  is  an  excellent  point  of 
departure  for  more  advanced  study  of  filtration,  and  that  their  efforts  are 
a  vast  improvement  to  the  totally  primitive  state-of-the-art  only  a  few  years 
ago.  It  is  equally  clear,  however,  that  significant  improvements  are  required 
if  filtration  is  to  be  put  to  the  best  possible  use  by  the  Navy. 

III.  DEVELOPMENT  OF  THE  THEORY 

This  section  outlines  the  conceptual  development  of  the  theory.  Much 
of  the  mathematics  is  presented  although  the  detailed  equations  are  reserved 
for  the  Appendix.  The  basic  concepts  follow  those  briefly  described  in 
Section  II,  2.1.  The  reader  may  skip  parts  2,  3,  and  5  without  loss  of  con¬ 
tinuity  in  the  discussion. 


1.  The  Governing  Parameters 


The  filtration  medium  is  assumed  to  have  a  porosity  e  (volume  void 


fraction),  to  be  of  length  L,  and  to  be  composed  of  "grains"  as  collectors 

of  radial  size  R.  At  this  point  both  the  contaminant  entrained  particle  and 

the  collector  are  idealized  as  spheres,  although  in  the  case  of  the  collector, 

cylinders  (fibers)  could  as  well  be  used.  The  particles  have  radii  a  and 

density  p.  The  lubricant  has  viscosity  p,  density  p  and  enters  the  filter  at 

velocity  V.  Particles  are  ultimately  captured  by  the  filter  grain  collector 

and  "stick"  due  to  van-der-Waals  forces  of  molecular  attraction.  This  force 

-13 

is  characterized  by  the  so-called  Hamaker  constant  H  a  10  dyne-cm,  where 
the  units  refer  to  force  times  the  distance  of  separation. 


12 


2.  The  Unit  Cell 

See  Figures  2  and  3  and  the  discussion  of  II,  2.1  for  a  description  of 
the  unit  cell. 

2.1  Size  of  the  Unit  Cell 

Following  the  model  of  Happel  [19],  the  unit  cell  is  of  such  a  size 
that  the  porosity  of  the  unit  cell  is  the  same  as  the  porosity  of  the  total 
filter.  Hius  it  can  be  shown  that 

where  Rg  is  the  unit  cell  radius,  see  Figure  3.  Hie  total  filter  efficiency 
is  found  from  the  unit  cell  efficiency.  Unit  cell  efficiency  is  the  ratio  of 
the  number  of  particles  which  ultimately  collide  with,  and  are  captured  by, 
the  collector  to  the  total  number  of  particles  which  enter  the  unit  cell. 

2.2  The  Fluid  Flow  Pattern  in  the  Unit  Cell 

The  following  velocity  field  follows  directly  from  solution  of  the 
Stokes  flow  equations  for  slow  flow  past  a  sphere,  as  presented  in  several 
texts,  e.g.,  Ref . [20] .  Spherical  coordinates  are  used  with  r  and  9  representing 
the  radial  and  angular  coordinates  shown  in  Figure  3.  The  flow  is  symmetric 
about  the  vertical  axis.  Hie  cp  angular  coordinate  rotates  around  the  vertical 
axis,  thus  the  cp-direct  ion  is  normal  to  the  plane  of  Figure  3,  not  shown. 

Symbols  v  and  vQ  denote  velocity  components  in  the  respective  directions: 

T 

,K  R3  K  R  K  r2 

v  -  -  Vcos  91 -=r-  +  +  K„  +  -2— 

r  "  r  ^  R2 

,  K  R3  K  R3  2K . r^\ 

v  -  V  .to  e(-  — r  ♦  -5-  ♦  v  — r-j  •  («> 

2r  R 

The  K’ s  are  functions  of  the  porosity  e,  with  the  exact  equations  presented  in 
the  Appendix. 
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2.3  Force  and  Torque  (Moment)  Balance  on  the  Particle 

Forces  and  moments  must  be  balanced  at  all  times  on  the  suspended 
particle  in  each  of  three  coordinate  directions.  Due  to  the  symmetry  of  the 
problem  forces  only  act  in  the  r-  and  6-directions  and  the  torque  vector  points 
only  in  the  ({^direction,  i.e. 


£F  *0  EFo»0  ZT  -0  .  (5) 

r  o  cp 

Likewise  all  particle  translation  occurs  in  the  r, 6-plane  and  the  rotation 
vector  points  in  the  cp-direction.  Forces  are  due  to  particle  inertia,  van- 
der-Waals  effect  (intermolecular  adhesion),  and  viscous  fluid  effects.  Torques 
sure  due  solely  to  viscous  effects.  By  nondimens ionalizing  the  equations  (5), 
as  performed  in  the  Appendix,  it  is  shown  that  the  magnitude  of  inertia  forces 
is  governed  by  a  Stokes  number 


St 


l  _£ 

9  p.R 


a2V 


=  2  i  f*') 

9  p  W 


Re 


(6) 


where  Re  is  the  Reynolds  number 


Re 


£VR 

V- 


(7) 


The  Reynolds  number  is  already  assumed  to  be  small  so  that  the  creeping  flow 
assumptions  can  be  used,  Section  III,  2.2.  Reynolds  number  is  as  a  general 
rule  negligible  in  lubrication  flows.  Hie  ratio  a/R  (particle  size  to  col¬ 
lector  size)  is  probably  less  than  one  although  the  particle  density  is  probably 
greater  than  the  fluid  density  (p  /p  >  1).  It  appears,  as  a  rule,  that  if  the 
Reynolds  number  is  small,  the  Stokes  number  is  also;  and  inertia  forces  can 
be  neglected. 
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2.4  The  Viscous  Fluid  Forces  and  Torques 

The  viscous  effects  are  extremely  complex.  This  is  due  primarily  to 
the  fact  that  while  simple  expressions  result  for  a  particle  in  an  unbounded 
fluid,  for  a  particle  near  an  obstacle,  many  complex  factors  come  into  play. 

For  example,  a  fluid  moving  at  velocity  u  in  a  motionless  unbounded  fluid  causes 
only  an  opposing  force  in  the  same  direction  by  "Stokes  law": 

F  =  6tt^  au  .  (8) 

If  an  obstacle  is  now  placed  nearby,  the  flow  field  around  the  particle  is  no 
longer  symmetric.  This  produces  unequal  pressures  and  shear  about  the  particle 
resulting  in  forces  in  both  directions  and  a  torque. 

In  the  absence  of  fluid  and  particle  inertia,  the  entire  problem  of 
evaluating  the  viscous  fluid  effects  is  linear  and  solutions  seprately  available 
in  the  literature  can  be  superimposed,  i.e.. 


The  vector  force  on  a 
translating,  rotating 
particle  in  a  shear¬ 
ing  fluid  flow  with 
vorticity  (fluid  ro¬ 
tation) 


The  vector  force  on  a 
translating  particle 
in  a  stationary 
fluid  [21] 


The  vector  force  on  a 
rotating  particle  in 
a  stationary  fluid  [22] 


The  vector  force  on  a 
+  stationary  particle  in 
a  shearing  fluid  [23] 


The  vector  force  on  a 
+  stationary  particle  in 
a  rotating  fluid  [23]  . 

~2  (9) 


■Hie  precise  technique  for  this  superposition  process  is  shown  in  the  text  by 
Rappel  and  Brenner  [20] . 


The  F^  and  F2  refer  to  vector  forces  on  a  moving  particle  in  a  stationary 
fluid.  The  radial  and  tangential  components  of  these  forces  (denoted  by  sub¬ 
scripts  r  and  @)  depend  on  the  radial  and  tangential  components  of  the  particle 
velocity  uy  and  Ug,  and  the  particle  rotation  o>: 
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F 


,  -»  F,  (u  ,uQ,ou) 

lr  lr  r’  9’ 


(10) 


pie  ■  ''le'W"1 

p2r  '  ^'W" 

P28  ’  • 

The  and  f 2  refer  to  vector  forces  on  a  stationary  particle  in  a  moving  fluid. 

The  radial  and  tangential  components  of  these  forces  depend  on  the  components 

of  fluid  velocity  v^  and  Vg,  at  the  particle  location  (rp,  9p).  The  fluid 

velocity  is  given  by  Eq.  (4),  now  evaluated  at  r=r  ,  9  =  9  : 

P  p 


fIr-£lr[Vr(Vep>’  V8  W  ‘  £lr<VV 


etc.  (11) 

The  eight  equations  forming  (10)  and  (11)  also  depend  on  parameters 

listed  in  Section  III,  1.  :  e,  R,  a,  p.,  and  V.  Densities  p  and  pp  no  longer 

have  influence  if  inertia  can  be  neglected.  The  van-der-Waals  constant  H 

does  not  affect  the  viscous  forces  and  torques.  The  total  filter  depth  L  enters 

later  in  the  total  efficiency.  For  the  unit  cell  the  depth  is  already  fixed 

by  ,ne  cell  radius  R  ,  Eq. (3).  The  exact  expressions  (10)  and  (11)  are  given 

s' 

in  the  Appendix.  Equation  (9)  now  has  the  form; 


Pvr  -  (ur,  ue,  w,  rp,  9p ; a,  e,  R, p,  V) 
Fv9  "  Fve ( W  ^  V  9p;  a?  e’ v) 


(12) 


The  semicolon  separates  the  parameters  from  the  independent  variables.  A 
similar  expression  is  also  derived  for  the  viscous  torques  T  on  the  particle, 
based  on  the  same  procedure  of  Eqs.  (9)  -  (12). 


y 


2.5  The  Force  and  Torque  Balance  Equations 


Equations  (5)  now  can  be  expressed  in  the  form: 
SPr“°“  Fvr(W<“» V0p?a’e>s^>v)  +  J?r(rp;R,H) 

Z  p0  “  0  “  Fv0  ( V  u0 >  ^  V  ®P '  a»  ® ’  R» ^  ’ V> 
2Tcp“0aTv9(WU)’V9p;a,e,R,lA,V)  ’ 


(13) 


where  3  is  the  known  van-der-Waals  force  which  only  occurs  in  the  radial 
r 


direction,  and  produces  no  torque. 


2.6  Solution  of  the  Force  and  Torque  Balance  Equations 

Hie  system  (13)  is  three  equations  in  three  unknowns  u^,  Ug,  and  u>  which 
can  be  solved  algebraically  as: 

ur-ur(rp,0p;a,e,R,n,v) 

u0“u0  (V  Va’e'R’tA'V) 


(14) 


<u  -  ou(r  ,0  ;a,e,R,n,V)  . 

p  p 


2 . 7  Nondimensional  Form 


Equations  (14)  can  be  written  more  concisely  as  a  nondimensional  set: 


I  -u  (r  ,0  ;  a,  e.  A) 
r  r  p*  p  ’  ’ 


WVV3,e,A) 


(15) 


S-‘"(rp,0p,a,e,A)  , 


where  r 

r  -  dimensionless  radial  particle  coordinate  ■ 
p  u  „  R 

u  Q  *  dimensionless  velocity  ■  — — 
r,  9  V 


tuR 


cu  *  dimensionless  particle  rotational  speed  »  — 


a  ■  dimensionless  particle  size  (relative  to  collector)  »  — 

A  -  Viscous  Effect  Number,  the  ratio  of  fluid  viscous  forces 
to  van-der-Waals  forces 


(16) 


9tt^R2V 


H 
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The  porosity  e  and  the  angle  0^  are  already  dimensionless.  A  dimensionless 
time  t  has  not  yet  been  used  but  will  be  required  shortly, 

t-H  . 


2.8  The  Particle  'trajectories 


The  coordinates  of  a  particle  r  and  0  vary  with  time  t  as  it  flows 

P  P 

along  its  trajectory.  The  particle  coordinates  are  related  to  the  velocities 


(18a) 


r*  if  • 


(18b) 


Dividing  (18b)  by  (18a)  and  rearranging  gives 

d6o  S0 

»  t a  g(r  ,0  ;a, e.A)  ,  (IS 

dr  r  u  p'  P  *  ’ 

Mr  r  4 

subject  to  some  initial  condition  0*0  when  r  »  r  . 

p  po  p  po 

The  solution  of  the  differential  equation  (19)  is  an  expression 

V  VV?po’V5,c,A)  (2C 

which  is  the  trajectory  of  a  given  particle,  see  Figure  3. 

Equations  (10)  through  (20)  are  extremely  complicated  and  the  opera¬ 
tions  ere  all  performed  by  computer. 

2.9  The  Limiting  Trajectory  and  Unit  Cell  Efficiency 


To  determine  the  efficiency  of  the  unit  cell,  it  is  not  necessary  to 

determine  the  trajectory  of  each  particle  which  enters  the  cell,  it  is  only 

necessary  to  find  the  trajectory  of  that  particle  which  is  just  barely  captured. 

The  particle  which  is  barely  captured  is  that  which  ends  up  at  0  •  rr  when 

po 


(21) 


F  >1.  see  Figure  3.  This  part..cle  enters  the  unit  cell  at  0  *9  and 
po  *  p  s 

-1/3 

r  »R/R-  (1-e)  ,  c.f.  Eq.  (3).  Substituting  in  (20)  gives 

p  s 

9s-0p[a-er1/3,  l,TT,I,e,A]  , 

where  0  is  the  angle  at  which  the  limiting  trajectory  crosses  the  cell 
s 

boundary. 

Any  particle  entering  the  cell  inside  (closer  to  the  vertical  axis 
than  this  "limiting  trajectory")  is  captured,  and  any  entering  outside,  escapes. 
It  is  assumed  that  particles  entering  the  cell  are  dispersed  uniformly  through¬ 
out  the  fluid.  Therefore  the  efficiency  1)  of  the  unit  cell  is 


Volume  flow  rate  of  fluid  entering  the  .  Q  .2 

. . .  .....  .  tt  (R  sin  o  )  V 

_  cell  xnside  the  limiting  traiectory _  =  s _ s 

o  ”  Total  volume  flow  rate  of  fluid  entering  rrp2v 

the  unit  cell  s 


(22) 


Thus  a  fairly  simple  expression  for  results: 

T)o  =  sin20s  =»  T^(a,  e,A)  .  (23) 

The  resulting  unit  cell  efficiencies  are  plotted  as  a  function  of 
the  relative  p article  size  a,  porosity  e,  and  the  viscous  effect  number  A  in 
Figure  4.  Physically  realistic  values  of  these  parameters  are  chosen  in  all 
cases,  as  discussed  below. 


3.  Total  Efficiency 

Total  efficiency  is  related  to  unit  cell  efficiency  in  a  straightforward 
way.  In  a  filter  of  length  L,  the  number  of  unit  cell  layers  N  is 

Filter  Depth  Number  N  »  ,  (24) 

s 

recalling  that  2Rg  is  the  diameter  of  the  unit  cell,  Eq.  (3),  and  hence  the 
thickness  of  one  layer.  If  n^(a)  particles  enter  the  filter  at  (relative) 


2 


size  a,  (a) [1  -  T|o (a)]  particles  escape  to  the  second  layer.  Of  these  particles 

(a)  [1  -  (a)]  [1  -  T)o  (a > ]  escape  to  the  third  layer,  etc.  After  N  layers,  the 

N 

number  of  particles  escaping  is  [1  -  T|q]  ,  and  the  number  captured  is 
n^  -  n^fl  -  Urns  the  total  efficiency  is 

Tl(S,e,A,N) -1-  [l-Tlo(a,e,A)]N  .  (25) 


4.  The  Frequency  Distribution  Function 

Although  a  number  of  systems  for  specifying  the  distribution  of  particle 
sizes  in  a  sample  of  particulate  matter  have  found  use  in  certain  fields,  the 
frequency  distribution  function  is  most  commonly  used  and  is  the  most  scienti¬ 
fically  valid.  It  can  be  visualized  in  the  following  way.  Consider  a  sample 
with  a  total  number  of  particles  nfc  in  some  radial  size  range,  say  5  to  50  pm. 
Assume  further  than  n^/n^. «  40%  of  the  particles  are  between  5  and  10  pm, 
n2/nt»30%  between  10  and  20  pm,  n-j/n^  =*  20%  between  20  and  30  (Jim,  and  n4/nfc »  10% 
between  30  and  50  pm.  The  size  increments  are  Aa^  »  5  pm,  Aa^  *  Aa^  =»  10  pm, 

Aa4  ■  20  pm.  The  midpoints  of  the  size  increments  are  =  7 . 5  pm,  a.^  =  15  pm, 
a3  *  25  pm,  a4  »  40  pm.  The  frequency  distribution  functions  are  defined  by 

An  ^  An2  ^ 

f(al)  f(a2)  ’  ’  StC-  (26) 

Assume  now  a  curve  is  fitted  through  a  plot  of  ffa^),  ...  ,  f(a4>  vs.  a  ,  a2, 
a3,  a4>  This  curve  is  the  frequency  distribution  function,  see  Figure  5. 

In  terms  of  calculus 

f<*>  TT’  <27> 

where  dn(a)  denotes  the  number  of  particles  in  a  small  size  increment  about 
size  a.  For  a  finite  size  increment 


f  (a)da 


(28) 


n  (a^  <  a  <  a2 )  “  nt  J 


For  the  entire  size  range  considered  a.  <  a  <  a. 

s  U 


and 


n(ag  <  a  <  «L>  • 


J  f(a)da«l 

a 


(29) 


For  all  possible  sizes  a  -  0  and  a  »  ® . 

S  u 


5.  the  Beta  Ratio 

From  Eqs. (27)  and  (28) 


n.(ao  <a  <»)  -  J  d^-n^  J  (30) 

a  a 

o  o 

where  fiD  (a)  is  the  particle  size  distribution  of  AC  Fine  Test  Dust,  at  the 
filter  inlet. 

From  the  definition  of  filtration  efficiency,  Eq. (1);  and  Eqs. (27) 
and  (28), 

dnfl-  tl~  T|(a)]dn^  ■  [1  -  TUa)]  n^f^  (a)da  , 


w 

V*o  <*  J  [i-noniDc 


)da 


(32) 


Therefore  the  beta  ratio,  defined  by  Eq. (2)  is 


0<a  ) 


n.  (a>a  )  n,  (a  <a  <®) 
i  o  i  o 


o  n  (a  >  a  )  n  (a  <  a  < »)  a> 
e  o  e  o  * 

J  El-  11(a)  ]fiD(a)da 


(33) 


Although  the  beta  ratio  has  been  defined  by  Fitch  et  al.;  in  words,  many  times 
in  various  articles,  apparently  this  is  the  first  rigorous  mathematical 
definition. 


IV.  RESULTS  OF  THE  FILTRATION  THEORY 

1.  Efficiency  Curves 

There  are  such  a  large  number  of  governing  parameters  for  total  efficiency, 
that  no  single  figure  can  adequately  portray  filtration  behavior.  Essentially 
one  figure  cam  only  show  the  effect  of  two  variables  on  efficiency  while  the 
others  are  held  constant.  The  horizontal  axis  in  all  figures  is  particle 
diametric  size  (major  dimension)  d  ■  2a,  in  micrometers.  In  the  analysis  of 
Section  III,  particle  size  relative  to  collector  size  D  is  used  (a-a/R*d/D) 
instead.  To  lend  physical  meaning  here,  a  value  for  collector  size  D  is 
assumed,  in  which  case  particle  size  d  has  physical  dimensions.  This  has  the 
net  effect  of  adding  a  dependent  parameter  to  the  efficiency  curves. 

The  governing  parameters  d  and  D  and  the  dimensionless  parameters 
porosity  6,  viscous  effect  number  A,  and  filter  length  number  N  were  defined, 
and  discussed  from  a  mathematical  standpoint  in  Sections  III,  1.,  Ill,  2.7, 
and  III, 3.  Their  significance  will  be  discussed  again  below  in  1.1  to  1.4. 

All  of  the  trends  below  will  be  seen  to  make  physical  sense. 

1.1  The  Effect  of  Filter  Length  N 

In  Figure  6,  the  effect  of  filter  length  on  the  efficiency  curves  is 

shown.  Porosity  is  held  constant  at  the  value  of  0.4  and  the  filter  grain 

10 

size  D  at  0.4  mm.  The  viscous  effect  number  A  is  held  constant  at  10  .  The 

filter  depth  parameter  is  N,  the  number  of  layers  of  collectors.  The  number 
of  layers  is  the  filter  length  divided  by  the  thickness  of  each  layer  (unit 
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cell),  the  mathematical  definition  for  N  is  Eq. (24).  It  seems  that  the  longer 
the  filter  the  greater  the  efficiency  and  the  steeper  the  efficiency  curve, 
which  seem a  physically  reasonable,  the  length  parameter  N  is  varied  from  10 
to  40. 


1.2  the  Effect  of  Porosity  e 

Porosity,  or  volume  void  fraction,  also  has  a  marked  effect  on  effi¬ 
ciency,  Figure  7.  As  the  porosity  goes  down,  filtration  efficiency  goes  up, 
and  the  lower  the  porosity  title  steeper  the  efficiency  curves.  As  the  filter 
becomes  more  and  more  densely  packed,  the  efficiency  (i.e.,  the  capture  per¬ 
centage)  goes  up,  which  again  makes  sense  physically. 

1.3  Fluid  Viscous  Effects  A 

the  parameter  A  is  defined  by  the  relation 

A  -  36  ,  <34, 

where  li  is  fluid  viscosity,  V  is  the  fluid  velocity  entering  the  filter  and 

-13 

H  is  the  van-der-Waals  adhesive  force  constant  as  10  ergs.  The  magnitude 
of  A  indicates  a  mean  ratio  of  fluid  viscous  forces  to  van-der-Waals'  molecular 
adhesive  forces  (between  the  particle  and  collector).  Note  from  Figure  8  that 
as  fluid  viscous  forces  go  up  (either  by  increased  viscosity  or  flow  velocity) 
filtration  efficiency  goes  down,  and  eventually  the  efficiency  curve  becomes 
very  broad  and  flat.  Since  van-der-Waals'  forces  tend  to  capture  particles 
it  seems  physically  reasonable  that  as  fluid  viscous  forces  overwhelm  the 

9 

capture  forces,  the  efficiency  will  go  down.  The  values  of  A  shown  (10  , 

10  11  12 

10  ,10  ,  and  10  )  are  physically  reasonable  for  actual  filters  and 


lubricants 


efficiency 
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Porhcle  Size  -  d  (/**») 

Figure  6  Theoretical  Efficiency  Curves  - 
The  Effect  of  Filter  Length 


Particle  Size  el  fyun) 


Figure  7  Theoretical  Efficiency  Curves  - 
The  Effect  of  Porosity 


Particle  Sue  d  (/***) 


Figure  8  Theoretical  Efficiency  Curves  -  Figure  9  Theoretical  Efficiency  Curves  - 
Fluid  Viscous  Effects  The  Effect  of  Filter  Grain  Size 
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1.4  The  Effect  of  Filter  Grain  Size  D 

As  the  filter  collector  or  grain  size  decreases,  the  efficiency  in¬ 
creases,  provided  the  other  parameters  remain  the  same.  Figure  9.  To  maintain 
constant  porosity,  the  grains  would  be  moved  closer  together  as  grain  size 
decreases.  This  effect  will  tend  to  trap  more  particles. 

2.  Beta  Ratio  Performance  Curves 

The  predicted  beta  ratio  performance  curves,  from  Eq. (33),  are  presented 
in  Figures  10-12.  The  overall  trends  are  very  much  like  those  for  the  effi¬ 
ciency  curves  and  will  not  be  discussed  in  detail.  Recall  that  the  beta  per¬ 
formance  curves  are  based  on  the  specific  particle  size  distribution  of  AC 
Fine  Test  Dust. 

these  figures  can  be  compared  to  Fitch's  Beta-ten  model,  refer  to 
Figure  13.  the  basic  form  of  the  theoretical  curves  here  and  the  Beta- ten 
empirical  curves  are  quite  similar.  Note,  however,  that  the  Beta-ten  model 
predicts  one  and  only  one  set  of  curves  for  all  filters.  None  of  the  effects 
listed  (1.1  to  1.4  above)  have  any  effect  whatsoever. 

V.  APPLICATION  OF  THE  THEORY  TO  A  LUBRICATION  SYSTEM 

This  section  outlines  an  application  of  the  filtration  theory,  just 
presented,  to  a  lubrication  system.  Two  simple  conclusions  are  drawn  which 
can  be  put  to  use  by  the  Navy  in  shipboard  systems  very  quickly.  The  conclu- 

1)  The  ideal  optimum  filter  is  " absolute"  with  the  cutoff  at  some  size 

d  (*  2a  ) .  Since  real  filters  are  not  absolute  the  best  filter 
o  o 

possesses  a  steep  efficiency  curve  at  some  size  dQ,  see  Figure  14. 

2)  The  optimum  size  cutoff  at  dQ,  for  a  given  lubrication  system, 
can  be  found  by  a  fairly  simple  series  of  field  tests. 


sions  are: 


1.  A  Simple  Conceptual  Model  of  a  Lubrication  System  with  Filtration 

Figure  15  depicts  schematically  a  simple  lubrication  system,  consisting 
of  an  engine,  pump  and  a  filter.  No  by-pass  is  shown  here  so  that  the  basic 
points  can  be  most  simply  illustrated,  although  this  and  other  features  can  be 
easily  added  to  the  model.  Particles  at  many  sizes  enter  at  the  engine,  either 
generated  by  on-going  wear,  or  due  to  ingression  of  dirt,  etc.  The  lubricant 
with  these  suspended  particles  is  pumped  around  to  the  filter.  Some  fraction, 
at  each  size  increment,  is  captured  by  the  filter  and  removed  from  the  flow, 
while  the  remainder  returns  to  the  engine  to  cause  additional  wear  and  other 
problems. 

2.  Analysis  of  the  Lubricant  System 
2.1  Particulate  Contamination  Content 

Total  particles  generated  in  the 
engine  by  wear  or  by  ingression 
(per  unit  volume  of  lubricant, 
per  unit  time) 

Particles  generated  in  the  engine 
and  leaving  at  size  a,  cf. 

Eq. (27). 

Number  of  particles  filtered  at 
size  a 

Number  of  particles  which  escape 
the  filter  and  return  to  the 
engine,  at  size  a. 


n 

(35a) 

-#» 

dn(a)  =  nf(a)da 

(35b) 

Tl(a)dn  (a)  *  n71(a)f  (a)da 

(35c) 

[1-  71(a)]dn(a)  * 

n  [1  -  T]  (a)]  f  (a)da  (35d) 


* 

As  before  dn(a)  denotes  the  number  of  particles  at  size  a  in  a  small  incre¬ 
ment  da. 


Vj 


I 


N 

'to, 

V- 
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Particles  captured* 


Figure  15  Schematic  of  a  Simple  Lubrication  System 
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2.2  Maintenance  Coats 


Now,  imagine  that  each  particle  carries  a  little  price  tag  Cf(a)  which 
represents  a  cost  to  filter  it.  Naturally  C^(a)  decreases  as  a  increases  - 
it  is  cheaper  to  filter  big  particles,  see  Figure  14.  Also  imagine  that  each 
particle  returning  to  the  engine  carries  another  price  tag  which  represents 
its  damage  potential  due  to  wear  C^fa).  The  Cw(a)  will  most  likely  increase 
as  a  increases  -  big  particles  cause  greater  wear,  see  Figure  14.  Therefore; 


Cost  of  filtering  particles  at 
size  a  (■  number  of  particles 
filtered  X  cost/particle) 

Cost  due  to  wear  of  particles 
at  size  a  (-  number  of  par¬ 
ticles  entering  engine  X 
cost/particle ) 

Total  cost  due  to  particles 
at  size  a  (sum  of  (a)  and 
(b)  above ) 

Total  costs  due  to  filtering, 
all  particle  sizes  [integral 
of  (a)] 

Total  costs  due  to  wear,  all 
particle  sizes  [integral  of 
(b)] 

Total  cost  of  lubrication 
system  maintenance,  for  all 
particle  sizes  [sum  of  (d) 
and  (e ) ] 

The  expression 


T](a)dn  (a)Cf  (a)  - 

nf  (a)T|(a)Cf  (a)da  (36a) 

[1-  Tl(a)]dn(a)C  (a)  - 
w 

nf  (a)Cw(a)  [1  -  Tl(a)]da  (36b) 

nf  (a){cf  (a)Tl(a)  + 

Cw(a)  [1-  11(a)]  ]da  (36c) 

CD 

F  -  J  nf  (a)T)(a)Cf  (a)da  (36d) 

o 

ao 

W  -  |  nf  (a)Cw(a)  [1- T1(a)]da  (36e) 

o 

00 

C  -  [  nf  (a)  {cf  (a)t](a)  + 

°  +  Cw(a)  [1- 11(a)]  }da  (36f) 


(36  f )  can  be  rearranged  to  read 

00  CD 

C«  n  J  f  (a)Cw(a)da  +  n  J  f  (a)T](a)  [Cf  (a)  -  Cy  (a)  Jda 


(37) 
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2.3  Minimizing  Maintenance  Costs  by  Optimizing  Filtration 

The  desire  is  now  to  minimize  maintenance  costs,  i.e.,  minimize  the 

expression  for  C,  Eg. (37),  by  selecting  the  best  filtration  performance  curve 

T)(a).  The  first  term  is  fixed  with  respect  to  the  filtration  and  cannot  be 

affected.  Thus  the  minimum  of  C  is  determined  solely  by  the  minimum  of  the 

second  term.  In  the  integrand  of  the  second  term,  f(a)  and  T)(a)  are  always 

positive  and  TJ (a)  is  between  zero  and  one.  The  bracketed  quantity  [Cf  (a)  - 

Cw(a)]  may  be  positive  or  negative,  see  Figure  14. 

Therefore  the  minimum  of  the  second  term  (and  consequently  the  minimum 

of  C)  is  obtained  by  maximizing  7] (a)  when  [Cf  (a)  -Cw(a)]  is  negative  and 

minimizing  T)  (a)  when  [Cf(a)  -Cw(a)]  is  positive.  Checking  Figure  14, 

[C _ (a)  -  C  (a)]  changes  sign  from  negative  to  positive  at  some  "cutoff’’  particle 
f  w 

size  aQ.  Therefore  the  optimal  filtration  performance  curve  if  (a)  is: 

if  (a)  »  0  a  <  ao 

7f  (a)  -  1  a  >  a  (38) 

o 

[C_  (a  )  -  C  (a  )]  -  0 
to  wo 

as  shown  on  the  figure. 

3.  Finding  the  Optimum  Filter  from  Lubrication  System  Tests 

From  the  discussion  of  Part  2,  it  appears  that  the  optimum  filter  is  one 

that  provides  a  cutoff  of  all  particles  greater  than  some  size  aQ.  The  size  aQ 

can  be  found  provided  the  functions  C  (a  )  and  C  (a  )  are  known.  This  part 

f  o  w  o 

outlines  a  simple  technique  for  finding  these  functions  from  lubrication  system 
tests,  perhaps  performed  on  board  ship. 

3.1  Required  Field  Data 

In  what  follows  it  is  assumed  that  a  means  to  systematically  vary  the 
level  of  filtration  exists,  for  example  by  inserting  different  filter  elements. 


In  addition  it  is  assumed  that  particle  size  frequency  distribution  functions 
in  the  lubricant  flow  can  be  obtained,  with  the  aid  of  some  particle  size 
analysis  device  or  "counter,”  and  that  filtration  efficiency  curves  are  available, 
either  by  laboratory  bench  tests  or  measured  directly  in  the  field.  The  latter 
procedure  can  also  be  readily  performed  with  a  particle  counter. 

Assume  that  this  data  has  been  taken  for  a  series  of  different  filters 
#1,  2,  . . .  ,  N.  Each  filter  causes  a  characteristic  particle  size  frequency 
distribution  function  to  occur  in  the  lubricant  flow  exiting  the  engine,  see 
curves  f  ,  f2,  ...  on  Figure  16.  As  a  rule,  the  finer  the  filter  the  higher 
the  f-curve  at  low  particle  sizes;  and  the  coarser  the  filter,  the  flatter  the 
frequency  curve.  Recall  that  the  curves  are  normalized  so  that  the  area  under 
the  curves  equals  one.  In  this  way,  frequency  distribution  function  charac¬ 
terizes  the  relative  distribution  of  the  particles  while  the  absolute  level  of 
contamination  is  determined  by  the  total  particle  count  n. 

For  each  filter,  certain  economic  data  has  been  taken: 

(1)  Total  costs  associated  with  using  each  filter  (cleaning,  filter 
replacement,  service  downtime  costs,  etc.).  These  are  the 
values  F  from  Eq. (36d)  -  F^,  F2,  ....  F^. 

(2)  Total  costs  associated  with  wear  of  the  entire  mechanical  system 
for  each  filter  (such  as  parts  replacement,  failure  downtime, 
etc.).  In  the  case  of  preventive  maintenance,  when  a  worn  com¬ 
ponent  is  changed  after,  say,  an  estimated  90%  of  its  useful  life, 
pro-rated  costs  should  be  estimated.  These  values  are  the  W  from 
Eq. (36e)  -  W, ,  W_  ...  w„.  The  values  F  and  W  are  shown  on 


Figure  16 


Frequency  Distribution  f(a) 
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Figure  16  Determining  Optimum  Filtration  from  Field  Data 
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3.2  Finding  the  Curves  Cf  and 

The  task  now  is:  given  F,  n,  f(a),  and  T|(a)  for  several  tests  1,  2 , 
...  N;  find  the  function  Cf  (a)  [ see  Eq. (36d)].  The  Cf  (a)  stays  the  same  for 
each  test.  The  same  problem  applies  to  finding  the  Cw  (a)  from  Eq. (36e). 
Taking  the  derivative  of  (36d) 


~  *  nf  (a)7)(a)Cf  (a) 


where  dF/da  can  be  expressed  as 


Rearranging  (39)  gives 


dF  _  dF  df 
da  df  da 


dp  df 
.  df  da 
f (a)  *  nf  (a)Tl(a) 


Each  of  the  quantities  on  the  right-hand  side  are  now  known.  Those  in 
the  denominator  are  directly  measured.  The  df/da  is  easily  obtained  from  the 


slope  on  Figure  16  at  each  a  and 


F  -  F 
2  *1 


df  f2(a)-f1(a) 

from  Figure  16.  The  procedure  is  then  repeated  over  the  spectrum  of  particle 
size  a.  Strictly,  only  two  different  filters  1  and  2,  fairly  close  in  proper¬ 
ties  are  required,  but  in  practice  a  series  of  curves  1,  2,  ...  N  will  provide 
increased  reliability.  Similarly  to  Eq. (41) 


Cw(a) 


dW  df 
df  d 

nf  (a)  [1-  Tl(a)] 


3.3  The  Optimum  Filter 


Hie  intersection  of  the  resulting  curves  (a)  and  Cf(a)  thus  provides 

the  optimum  filter  size  a  .  Since  no  filter  is  absolute  at  size  a  ,  the  actual 

o  O' 
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test  filter  is  one  with  a  steep  slope  in  the  efficiency  curve  at  size  aQ, 

Figure  14. 

4.  Osing  the  Filtration  Theory  to  Design  the  Best  Filter 

Recalling  the  discussion  of  Section  IV,  1,  and  Figures  6-9,  a  sort  of 
"tuning"  can  be  achieved  by  adjusting  the  filtration  parameters  to  approach  the 
optimum  curve.  Note  that  increasing  filter  depth,  decreasing  filter  grain  size 
and  decreasing  porosity  all  serve  to  significantly  steepen  the  efficiency  curves 
(especially  the  latter  two).  The  "cutoff"  frequency  is  also  decreased  by  these 
same  parameters. 

therefore,  by  juggling  these  three  parameters,  a  filter  with  almost  any 
desired  efficiency  curve  could  be  constructed.  For  example,  assume  the  desired 
cutoff  is  20  nm.  An  iteration  process,  using  repeated  computer  runs  of  the 
filter  theory,  may  be  used  in  the  following  ways 

1)  Start  with  N»25,  A  =*10,  D*.4mm,  e»0.4;  the  middle  curves  on 

Figures  8  and  9,  this  curve  goes  through  20  |j,m  at  about  56% 
efficiency. 

2)  Steepen  the  curve  by  adding  greater  depth,  say  N»60.  this  effect 
is  shown  on  Figure  6,  but  the  specific  curve  for  N  **  60  is  not 
shown. 

3)  Step  2)  will  also  shift  the  curve  to  the  left,  so  now  it  may  go 
through,  say,  20  (im  at  70%  efficiency. 

4)  Shift  the  curve  back  to  the  right  by  decreasing  the  porosity, 
say,  to  0.35,  as  suggested  in  Figure  9; 

5)  Repeat  2)  to  4),  etc. 

Eventually  the  desired  filter  curve  can  be  approached. 
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Even  if  the  theoretical  predictions  are  not  entirely  accurate,  the  trends 
suggested  may  greatly  aid  designers  in  constructing  the  optimum  filter  for  a 
particular  application. 

VI.  EXPERIMENTAL  RESEARCH 

The  experimental  research  conducted  to  date  consists  of  measuring  the 
rheological  properties  of  used  oils  and  other  fluid/particle  systems,  obtained 
from  various  sources  under  various  conditions  of  service,  and  correlating  these 
properties  with  the  particulate  content  of  the  oil. 

1.  The  Fluids  Tested 

Five  fluids  have  been  tested  to  date.  In  all  cases  the  clean  oil  is  an 
identical  sample  of  comnercial  SAE-30  detergent  motor  oil.  The  commercial 
grade  oil  was  used,  rather  than  a  standard  base  mineral  oil  with  clearly 
specified  properties  and  content,  due  to  the  presence  of  the  detergent.  The 
detergent  serves  to  suspend  the  smaller  particles,  and  particulate- loaded  oils 
of  this  type  probably  better  simulate  an  actual  field  environment. 

The  five  fluid  samples  were: 

Fluid  1:  An  AC  Fine  Test  Dust  Suspension,  0.5  gr  particles/50  cm3 
fluid. 

Fluid  2:  A  used  oil  from  1000  miles  gasoline  engine  automotive 
road  service. 

Fluid  3:  A  used  oil  from  1000  miles  diesel  automotive  road  service. 

Fluid  4:  A  used  oil  from  30  hours  of  laboratory  dynamometer 
gasoline  engine  operation. 

Fluid  5:  A  used  oil  from  30  hours  of  laboratory  dynamometer 
diesel  engine  operation. 
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Oils  from  automotive  use  were  used  because  they  were  readily  available 
for  this  project,  and  because  some  degree  of  control  could  be  exerted  as  to 
type  of  oil,  length  and  type  of  service  and  engine  type.  Hie  1000  miles  of 
field  service  is  roughly  equivalent,  in  time,  to  30  laboratory  hours. 

2.  Determination  of  the  Particulate  Content 

Hie  particle  size  analysis  was  performed  using  the  H1AC  Model  PC  320 
Particle  Counting  and  Size  Analysis  System,  purchased  under  the  provisions  of 
this  grant. 

For  the  fluids  tested,  particle  counts  were  run  in  (diametrical)  size 
increments  of  2.5  pm,  from  2.5  tun  to  100  tun.  It  was  found,  when  preparing 
particle  size  frequency  distribution  functions,  that  the  great  preponderance 
of  particles  in  the  range  2,5-10  pm  tends  to  mask  the  important  characteris¬ 
tics  in  the  more  crucial  range  of  10  <  d  <  100  pm.  For  this  reason,  this 
lower  range  was  omitted  when  preparing  the  frequency  distributions. 

Hie  total  particle  counts  in  the  ranges  10<d<100  pm  and  2.5  <d  <100  pm 
are  presented  in  Table  I. 

TABLE  I 

TOTAL  PARTICLE  COONTS 

cm 


Fluid 

n  (10<d<100pm) 

11^2.5  <d  <  100pm) 

% 

n 

a 

1  (AC  dust) 

10,640 

72,819 

6.84 

2  (Gas,  field) 

105,525 

917,400 

8.69 

3  (Diesel,  field) 

52,170 

891,433 

17.1 

4  (Gas,  lab) 

25,533 

24,945 

9.77 

5  (Diesel,  lab) 

46,410 

589,407 

12.7 
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Note  that:  (1)  the  diesel  engine  oils  are  dirtier  than  the  gasoline  engine 
oils,  (2)  the  field  service  engine  oils  sure  dirtier  than  the  laboratory  engine 
oils,  (3)  diesel  engine  oils  have  a  much  higher  percentage  of  very  small 
suspended  particles,  and  (4)  AC  Fine  Test  Dust  has  a  much  lower  percentage  of 
very  small  suspended  particles. 

The  frequency  distribution  functions  for  the  five  samples  (10  <d  <100 (Am) 
jure  shown  in  Figure  17.  Note  that  although  all  curves  show  the  same  downward 
trends,  the  following  substantial  differences  jure  noted:  (1)  at  a  given 
particle  size,  one  oil  may  have  double  or  triple  the  particulate  content  of 
another,  (2)  the  content  of  AC  Fine  Test  Dust  is  substantially  different  from 
that  of  either  oil.  These  facts  are  illustrated  in  the  blow-up  of  the  curves 
in  the  20  to  40  (im  range,  also  shown  in  Figure  17. 

3.  Rheological  Measurements 

3.1  The  Capillary  Viscometer 

A  special  capillary  viscometer  instrument  has  been  constructed. 

Figure  18.  The  instrument  is  capj&le  of  measuring  viscosity  at  a  combination 
of  temperatures  and  shear  rate.  The  apparatus  consists  of  (a)  a  nitrogen 
pressure  supply  tank,  (b)  a  plenum  bottle  and  pressure  gage,  (c)  plumbing 
to  allow  addition  of  the  fluid  seunple  and  application  of  the  known  pressure 
in  the  plenum  tank  to  the  sample,  (d)  a  constant  temperature  bath,  (e)  the 
capillary  itself,  (f)  a  flow  measuring  device  which  electronically  times  the 
passage  of  a  soap  bubble  between  two  marks  which  enclose  a  known  volume,  and 
(g)  plumbing  to  push  the  fluid  sample  back  through  the  capilletry  to  perform 
another  run. 

ttie  viscosity  is  proportional  to  the  product  of  the  plenum  pressure 
and  the  bubble  passage  time,  subject  to  several  correction  factors  [24] . 

The  apparatus  is  calibrated  with  a  standardized  fluid  of  known  viscosity. 
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3.2  Measured  Results 

Curiously  the  viscosities  of  the  fluid  samples  (1)  to  (5)  were  essenti¬ 
ally  Newtonian,  i.e.,  viscosity  did  not  change  with  rate  of  shear.  This 
behavior  is  at  odds  with  a  great  deal  of  rheological  research  on  suspensions 
of  clay  particles  in  water,  slurries,  etc.  [25,26].  The  overall  viscosity 
level,  relative  to  that  of  the  base  oil  is  significantly  increased  due  probably 
to  chemical  contamination  and  sub-micron  size  suspended  carbonaceous  matter. 

The  larger  suspended  particles  appear  to  have  little  influence. 

The  viscosity  curves  are  shown  in  Figure  19.  Note  the  significant 
increase  in  the  used  oil  viscosity  levels  and  the  negligible  change  due  to 
the  AC  Fine  Test  Oust,  which  would  tend  to  indicate  that  the  viscosity  level 
change  is  not  due  primarily  to  the  particulate  level  itself. 

Perhaps  at  higher  shear  levels  greater  non-Newtonian  effects  may  appear. 
These  higher  levels  will  be  obtainable  with  a  new  capillary  series  now  being 
constructed. 

3.3  Continuing  Work 

Other  capillaries  are  being  constructed  which  will  allow  extension  of 
the  shear  rate  to  5  x  IQ5  sec  \  which  approaches  that  occurring  in  rolling 
contact  bearings.  In  addition,  the  larger  particles,  say  d  >  5  will  be 
filtered  from  each  sample  and  viscosity  curves  taken.  If  the  curves  are  still 
high  relative  to  the  base  oil,  the  conclusion  will  surely  be  that  viscosity 
change  is  due  to  minute  particles  and  chemical  contamination,  as  opposed  to 
the  relatively  larger  filterable  particles. 

Studies  will  continue  on  additional  fluid  samples  from  other  appli¬ 
cations.  Viscosity  data  and  particle  size  analysis  have  been  performed  on 
polymer-additive  multigrade  oils,  which  clearly  demonstrate  non-Newtonian 


Viscosity  (cp.) 
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flut'd  (3)  Diesel  flrf,  Field  b 
Fhtd  (4)  Gas  iat>  4~ 
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Figure  19  theological  Properties  of  Used  Oil 
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behavior  in  the  clean  and  contaminated  conditions.  The  significance  of  these 
tests  are  not  now  conclusive. 

4.  Conclusions 

important  experimental  conclusions  are  as  follows: 

1.  The  frequency  distribution  functions  of  used  oils  are  widely 
different  and  cannot  be  simply  characterized  by  that  of  AC  Fine 
Test  Dust. 

2.  The  frequency  distribution  functions  and  the  level  of  contamination 
of  filterable  particles  (5  urn  <  d)  has  little  influence  on  lubri¬ 
cant  flow  behavior,  at  least  for  nonmultigrade  oils.  This  finding 
may  have  implications  regarding  extending  useful  lubricant  life. 

VII.  FUTURE  RESEARCH 

As  discussed  in  the  Introduction,  the  originally  proposed  two-year 
research  program  consisted  of  four  tasks,  although  only  the  first  two  task 
areas  were  addressed  in  the  funded  one-year  effort.  The  remaining  two  areas 
are: 

(3)  Perform  controlled  experiments  on  filtration  to  test  the 
theory. 

(4)  Perform  centrifugal  separation  studies  on  used  oil,  with  a 
unit  currently  available  at  RPI,  in  conjunction  with  on-going 
wear  studies. 

In  view  of  the  success  of  the  filtration  theory  to  date  (Sections  III 
and  TV),  the  promise  of  a  quick  pay-back  to  the  Navy  in  applications 
(Section  V),  and  the  limited  resources  available,  the  principal  investigator 
feels  that  additional  studies  in  filtration  will  currently  be  more  profitable 
than  Task  (4),  the  centrifugal  studies. 

It  is  proposed  therefore  that  Task  (4)  above  be  omitted  and  that 
future  research  consist  of  two  task  areas: 
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(1)  Continue  the  filtration  theory  work.  Continuing  research 
will  be  broken  down  into  five  areas: 

(a)  Development  of  a  simpler  user  -  accessible  computer 
program. 

(b)  Study  of  deposit  build-up  and  clogging. 

(c)  Study  the  capture  of  flat,  flake-like  particles. 

(d)  Pressure  drop,  and  particle  inertia  studies. 

(e)  Study  of  fiber-like  collectors. 

(2)  Filtration  experiments  to  test  the  theory  described  above. 


1.  Additional  Work  on  the  First  Year's  Program 

In  the  remaining  two  months  on  the  first-year  program,  several  on-going 
activities  will  be  completed. 

1)  Develop  an  approximate  formula  for  filtration  efficiency.  A 
semi-empirical  equation  for  efficiency  T|  as  a  function  of  the 
relative  particle  size  a-,  porosity  e,  depth  N,  and  viscosity 
number  A  will  be  developed.  This  can  be  used  by  designers, 
etc. ,  although  with  reduced  levels  of  accuracy. 

2)  Continue  rheological  measurements  and  particle  counts  on  used 
oil  samples.  This  item  was  discussed  under  VI,  3.3. 

2.  Future  Research  Program  -  Continue  Filtration  Theory  Studies 
2. 1  Development  of  a  Simple  User-Accessible  Program 

The  present  program  is  very  complicated,  expensive  and  time  consuming, 
requiring  a  large  computer.  The  efficiency  of  this  program  can  be  increased 
considerably.  A  simple  FORTRAN  program  listing  will  be  made  available, 
trimmed  to  100  lines  or  so,  that  can  be  put  on  a  small  computer  and  used  by 
lubrication  system  or  filtration  designers,  etc.  Eventually  it  is  hoped  that 


simplified  program  can  be  put  on  a  hand  calculator  card  insert. 
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2.2  Study  of  Deposit  Build-up  and  Clogging 

This  important  effect  can  be  included  in  an  improved  theory.  The 
deposited  particles  will  now  become  collectors  themselves.  Some  may  build  up 
within  the  filter  medium  decreasing  the  porosity  and  increasing  effective 
collector  size.  Eventually  a  dense  "cake"  may  build  up  on  the  filter  inlet 
side,  which  itself  becomes  a  filter  medium.  In  this  case  the  collectors  of 
the  cake  will  be  of  distributed  size  range.  This  phenomenon  is  very  important 
but  will  add  considerably  to  the  complexity  of  the  theory. 

2.3  Study  the  Capture  of  Flat  Flake-Like  Particles 

The  present  studies,  in  principle,  apply  to  particles  that  sure 
spherical,  or  nearly  so.  A  general  theory  for  arbitrary  shaped  particles 
is  not  feasible  at  this  time,  but  flat,  two-dimensional  flake- like  particles 
can  be  accommodated.  Between  these  two  extremes,  spheres  and  flakes,  the 
behavior  of  most  types  of  particulate  matter  in  the  filter  can  be  predicted. 
Fortunately  many  wear  particles  are  indeed  shaped  like  flat  flakes. 

2.4  Pressure  Drop  and  Particle  Inertia  Studies 

Simple  existing  models  for  pressure  losses  through  porous  media  will 
be  included  in  the  theory.  Certain  experimental  and  analytical  studies  to 
date  suggest  that  pressure  drop  predictions  are  unreliable  but  they  can  be 
checked  against  the  experimentally  obtained  values  (see  below).  Perhaps  the 
simple  models  can  be  modified  in  a  semi-empirical  fashion  to  adequately 
portray  pressure  drop  behavior.  Particle  inertia  forces  can  be  included  in 
the  momentum  balances,  cf.  Eq. (A. 8)  of  the  Appendix. 
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2.5  Study  of  Fiber- Like,  as  Opposed  to  Grain-Like,  Collectors 

This  part  is  a  very  straightforward  reworking  of  the  filtration 
theory  development.  The  fluid  velocities,  Eq. (4),  are  now  different  but  the 
entire  theory  proceeds  as  before,  from  Eqs. (5)  through  (25). 

3.  Future  Research  Program  -  Filtration  Experiments 

3.1  Experiments  to  Reproduce  the  theoretical  Setting 

Carefully  controlled  experiments  will  be  performed  using  spherical 
collectors  with  uniform  spacing,  once- through  flow,  and  spherical  suspended 
particles.  Particle  counts  will  be  obtained  up  and  downstream  of  the  filter. 
All  of  the  theoretical  parameters  will  be  precisely  measured.  The  emphasis 
will  be  on  reproducing  the  theoretical  conditions,  to  test  the  theory  and 
modify  the  theory  accordingly  if  experimental  data  so  suggest. 

3.2  Lubrication  Filtration  Experiments 

The  emphasis  here  is  to  create  a  more  real-life  set  of  conditions. 

Hie  flow  will  be  multipass.  Both  oil  and  AC  Test  Dust,  or  actual  used  oil, 
will  be  used  as  the  fluid/particle  system.  The  governing  parameters  will  be 
measured  as  accurately  as  possible  although  there  will  be  many  deviations 
from  the  theoretical  conditions,  e.g.,  collectors  will  be  nonspherical  of  a 
various  size.  Therefore  this  set  of  experiments  will  not  rigorously  test 
the  theory  itself  but  rather  indicate  whether  or  not  idealized  theoretical 
assumptions  can  be  reasonably  applied  to  a  real-life  filter  situation. 


46 


REFERENCES 


1.  Peterson,  M.B.  et  al.,  "Feasibility  Study  for  a  Diesel  Engine  Condition 
Honitoring  System  for  1179  Class  LSTs,"  Mechanical  Failures  Prevention 
Group,  Report  NBS-GCR  74-43,  August  1975. 

2.  Peterson,  M.B.  et  al.,  "Wear  Control  in  Naval  Aircraft,"  Naval  Air  Develop¬ 
ment  Center,  Project  Report,  ARP/SLP  Project  Office,  Warminster,  PA,  1975. 

3.  Sin,  H. ,  Saka,  N.  and  Suh,  N.P.,  "Abrasive  Wear  Mechanisms  and  the  Grit 
Size  Effect,"  Wear.  55.  163-190  (1975). 

4.  Avient,  B.w. E.  and  Wilman,  H. ,  "New  Features  of  the  Abrasion  Process  Shown 
by  Soft  Metals;  the  Nature  of  Mechanical  Polishing,"  Br.  Jour.  Appl. 
Physics,  13.  521-526  (1962). 

5.  Rabinowicz,  E.  and  Mutis,  A.,  "Effect  of  Abrasive  Particle  Size  on  Wear," 
Wear.  £,  381-390  (1965). 

6.  Larsen-Basse,  J. ,  "Influence  of  Grit  Diameter  and  Specimen  Size  on  Wear 
During  Sliding  Abrasion,"  Wear.  12,  35-53  (1968). 

7.  Davies,  C.N. ,  Air  Filtration,  Academic  Press,  London  1973. 

8.  Theodore,  L.  and  Buonicore,  A.J. ,  Industrial  Air  Pollution  Control  for 
Particulates.  CRC  Press,  Cleveland,  1976. 

9.  Payatokes,  A.C.,  Rajagopolen,  R.  and  Tien,  C.,  "Application  of  Porous 
Media  Models  to  the  Study  of  Deep  Bed  Filtration,"  Canad.  Jour.  Chem.  Eng., 
52,  722-731  (1974). 

10.  Rajagopalen,  R.  and  Tien,  C. ,  "Trajectory  Analysis  of  Deep-Bed  Filtration 
with  the  Sphere  in  Cell  Porous  Media  Model,"  A. I.Ch.E.  Jour..  22,  3, 

523-533  (1976). 

11.  Payatakes,  A.C.,  Tien,  C.  and  Turian,  R.M. ,  "Trajectory  Calculation  of 
Particle  Deposition  in  Deep-Bed  Filtration,"  A. I.Ch.E.  Jour..  20.  889-900 
(1974). 

12.  Loewenthal,  S.H.  and  Moyer,  P.W. ,  "Filtration  Effects  on  Bull  Bearing  Life 
and  condition  in  a  Contaminated  Lubricant. "  Jour,  Lub.  Tech.,  ASMS  Trans., 
101.  2,  171-180  (1979). 

13.  Fitzsimmons,  B.  and  Clevenger,  H.D.,  "Contaminated  Lubricants  and  Tapered 
Roller  Bearing  Wear,"  ASLE  Trans. .  20,  2,  97-107  (1977). 

14.  Lynch,  C.W.  and  Cooper,  R.  B. ,  "The  Development  of  a  Three- Micron  Absolute 
Main  Oil  Filter  for  the  T53  Gas  Turbine,"  Jour.  Lub.  Tech. .  AS  ME  Trans. . 

93,  3,  430-436  (1971). 


v 


47 


15.  Schilling,  A.,  Automobile  Engine  Lubrication.  Scientific  Publications 
(G.B. )  Ltd.,  Broseley,  England,  1972,  p.44. 

16.  Bensch,  L.E. ,  Fitch,  E.C.  and  Tessman,  R.K.,  Contamination  Control  for 
the  Fluid  Power  Industry.  Pacific  Scientific  Co.,  1978. 

17.  Fitch,  E.C. ,  "The  Multipass  Filter  Test  -  Now  a  Viable  Tool,"  Basic  Fluid 
Power  Research  Program.  S3,  Paper  P74-39,  Oklahoma  State  University,  1974. 

18.  Fitch,  E.C. ,  "Controlling  Contaminant  Wear  Through  Filtration,"  Wear,  34, 

3,  319-330  (1975). 

19.  Happel,  J. ,  "Viscous  Flow  in  Multiparticle  Systems:  Slow  Motion  of  Fluids 
Relative  to  Beds  of  Spherical  Particles,"  A. I.Ch.E,  Jour.,  4 ,  147  (1958). 

20.  Happel,  J.  and  Brenner,  H. ,  Low  Reynolds  Number  Hydrodynamics,  Prentice 
Hall,  Englewood  Cliffs,  New  Jersey,  1965. 

21.  Brenner,  H. ,  "The  Slow  Motion  of  a  Sphere  Through  a  Viscous  Fluid  Towards 
a  Plane  Surface,"  Chem.  Eng.  Sci. .  16.  242-251  (1961), 

22.  Goldman,  A.J. ,  Cox,  R.G.  and  Brenner,  H. ,  "Slow  viscous  Motion  of  a  Sphere 

Parallel  to  a  Plane  Wall  -  I.  Motion  Through  a  Quiescent  Fluid,"  Chem.  Eng. 
Sci..  22,  637-651  (1967).  " 

23.  Goren,  S.  and  0*Neill,  M.,  "On  the  Hydrodynamic  Resistance  to  a  Particle 
of  a  Dilute  Suspension  when  in  the  Neighborhood  of  a  Large  Obstacle," 

Chem.  Eng.  Sci..  26.  325-338  (1971). 

24.  Fredrickson,  A.G. ,  Principles  and  Applications  of  Rheology.  Prentice-Hall, 
Englewood  Cliffs,  New  Jersey,  1964,  pp. 16-23. 

25.  Metzner,  A. B.  and  Reed,  J.C. ,  "Flow  of  Non-Newtonian  Fluids  -  Correlation 
of  the  Laminar,  Transition  and  Turbulent  Flow  Regions,”  A. I.Ch.E.  Jour.. 

1,  434-440  (1955).  ' 

26.  Orr,  C.,  Particulate  Technology.  MacMillan  Company,  New  York,  1966, 
pp, 126-128. 


APPENDIX 


Al.  Purpose 

The  purpose  of  this  appendix  is  to  present  the  detailed  equations  and  steps 
of  the  theoretical  development,  which  were  put  forth  in  general  functional  form 
in  Section  III,  2.  The  equations  follow  with  relatively  little  explanation  as 
the  overall  concepts  were  discussed  adequately  in  the  main  text.  The  breakdown 
by  numerical  heading  parallels  that  of  Section  III. 


A2.  Governing  Equations 

A2.1  Size  of  the  unit  Cell 

Adequately  covered  in  Section  III,  2.1. 

A2.2  Fluid  Flow  Pattern 

The  K-parameters  follow: 


p  -  (1-e) 
W  »  2  -  3p  +  3p 


1/3 

5 


-1 


K2-PO+2P  )/W 

A2.2a  Other  Dimensionless  Variables 


K2=  (3  +p5)/W 
X4  -  -  p5/w 


z  -  ^-5  x-r(0-  0  ) 

a  p 


(A.  1) 


(A.  2) 


The  velocity  field  can  then  be  approximated  to  order  2  in  x  and  z  by 


-3—2 

v  ■  -  Ba  z  cos  0 
z  p 

_  _  _2_  _3_2 

v  ■  Ba  x  z  cos  0  +Ca  z  sin  0  +Di  z  sin  0  (A.  3) 

x  p  p  p 


where 


B(e) 


5 


3-3p~ 

W 


D(e) 


C  (e ) 

7  ,  5 

-  2  '  3P 

W 


3- 


w 


(A.  4) 


A2.3  Force  and  Torque  Balance  Equations 


1  o  t  r 
£  T  *0»t  +  tF  +  T  +T 

cp  <p  cp  cp  cp 

du 

r  s  n 

ZF  =  m  «F  +  f  +F 
r  dt  r  r  r 

svir-'e  fetF-8+p9 


(A.  5) 


A2.4  Viscous  Fluid  Forces  and  Torques 

The  lower  case  f  and  t  indicate  viscous  forces  (torques)  on  a  stationary 
particle  in  a  moving  fluid,  and  capital  F  and  T  indicate  viscous  forces 
(torques)  on  a  moving  particle  in  a  stationary  fluid.  The  code  for  the  super¬ 
scripts  is  as  follows: 

t  -  particle  translates  tangential  to  surface 
n  -  particle  translates  normal  to  surface 
r  -  particle  rotates  with  axis  parallel  to  surface 
Jt  -  linear  shearing  fluid  flow 
p  -  parabolic  shearing  fluid  flow 

s  -  axisymmetric  stagnation  fluid  flow  towards  the  collector. 

Xn  the  most  general  case, all  of  the  above  fluid  and  particle  motions  are  super¬ 
imposed.  The  fluid  particle  mass  is 


8  3 

*  "  3  pp  * 


(A.  6) 


The  adhesive  force  is  F?.  The  T  and  F  are  functions  of  the  fluid  velocities 
at  the  particle  center  (if  the  particle  were  not  there),  and  the  t  and  f  are 
functions  of  the  particle  velocities 
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t*«  t* 
cp  tp 


.1  3/VCS  sin  9v 

:cp8^a3( - 5— £) 


%  8"lia3(VDa2zpsin9p) 


T^aT^  8V3(») 

KmK  6^a  (VBaV  cos  9  ) 
r  r  P  p 

A  —A 

*  f  6rtvia  (VCaz  sin  9  > 

«  9  p  p 

f 0  *  fg  6TT\ia(VD522p  sin  9p) 
P^-pP  STTp.aCUg) 


f”*p£  6ttjj,  a(ur) 
Pg  *  Pg  6tt|x  a(o)V) 

P9**9  6ttU-  ®(uQ) 


(A.  7) 


r  3a (5  ~1)2[Z  +1)2 
P  P 

The  system  of  equations  (A. 7)  is  set  up  so  that  the  barred  symbols  are 

correction  factors  to  Stokes  law  for  forces  or  torques  on  a  particle,  due  to 

the  presence  of  the  obstacle  (collector),  in  an  unbounded  fluid,  or  as  z  -» 

’  P  * 

the  barred  symbols  approach  one.  These  correction  factors  are  often  quite 
significant,  for  instance  when  Sp-*1,  F-«.  The  quantity  in  parentheses  is 
the  pertinent  rotational  speed  (in  the  case  of  torque)  or  translational  speed 
which  appears  in  the  basic  form  of  Stokes  law,  Eq.  (18).  The  values  of  the 
barred  quantities  Tp,  ...  ,  fj  are  found  in  Refs. [10]  and  [21-23],  as  functions 
of  the  particle  radial  location  z  . 


Substituting  (A.  6)  and  (A.  7}  into  {A.  5)  and  re  air  ranging  to-  nondimen- 


sional  form  gives: 


%  tfsy  sin  0  +  Dt^i2z_  sin  0  +  uQTfc  +  off1  = 0 
‘  tp  P  P  tpp  p  9  cp  cp 


(z-ir  (z+i) 
p  p 


2  2  uu 

■s  +  Bf  a  z  cos  0  +  u  P-St  — - 
2  r  P  p  r  r  dt 


i  «  •  _  du 

cf9SIp  *in  V”1!*  5p  3i”  VVi**5!-*  IT 


A2.6  Solution  of  the  Force  and  Torque  Equations 


(A.  8) 


Equations  (A. 8)  can  be  solved  for  u  and  u  as 

r  0 

3  ■  |”a - - - -  +  Bf3®2^2  cos  0  1 

r  FL  (Z  -1)2(2  +1)2  r  P  P-l 

r  p  p 

ff9  '  93  >ln  6p'-C?H  '  "1*^  +f  V  +D52S\! 


(A.  9) 


(d)  is  not  required)  Where 
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(A.  10) 


A2.7  Mondlmensional  Form 

Recall  that  B,  C  and  D  are  functions  of  e  and  z  =*  a  (r  -  1). 

P  P 

A2.8  The  Particle  Trajectories 


y  m  g3^  8in  9P  r^V  (TP"1)~D^5  (?p~1)  +  2  v+py  <V1}  i 

drp  5o  LA  —3 - ^—3 -  +  Bf®34  (?  -  l)2  cos  0  -* 

tr®  -i>-i][r(2  -D+ij  r  p  p 

p  p 

(A.  11) 


A2.9  The  Limiting  Trajectox 


Initial  conditions  for  (A. 11)  are  0  *  tt  at  ?  » 1.  a  fourth-order 

P  P 

Runga-Kutta  scheme  is  used  to  solve  Eq. (A. 11)  with  the  incremental  values 


A5. 


of  Arp  determined  by  the  nature  of  the  local  trajectory.  Near  the  regions  0 »  tt 

and  9*0,  fairly  large  values  will  suffice,  Ar^  «  Rg/(100  R);  as  large  Ar^ 

produces  fairly  small  A6  ,  see  Figure  3.  Near  9  ■  tt/2,  AF  «  R  /  (10, 000  R). 

P  P  P  s 

The  correction  factors  F,  f,  T,  t  and  g  are  recomputed  at  each  step.  The 

procedure  continues  until  r  *  R  /R,  whence  6*9. 

p  s  p  s 


